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NONUNIQUENESS OE SOLUTIONS OE THE NAVIER-STOKES 
EQUATIONS ON NEGATIVELY CURVED RIEMANNIAN 

MANIFOLDS 

LEANDRO A. LICHTENFELZ 


Abstract. In a well-known work, M. Anderson constructed a Hadamard man¬ 
ifold which carries non-zero harmonic p-forms when p ^ nl2^ thus 

disproving the Dodziuk-Singer conjecture. In this paper, we use the manifold 
(M^^g) in order to solve another problem in geometric analysis, namely the 
nonuniqueness of solutions of Leray-Hopf type of the Navier-Stokes equations. 


1. Introduction 

The aim of this paper is to prove that the Cauchy problem for the Navier-Stokes 
equations is ill-posed on the Riemannian manifold {M^^g) constructed by Anderson 
[T] in his counterexample to the Dodziuk-Singer conjecture (see also [TD], [^, [5], 
M)- This answers a question raised by Chan and Czubak in [B], pursued also by 
Khesin and Misiolek |18) . 

In order to provide the appropriate context, we begin with a brief review of some 
pertinent results connected to the uniqueness problem. Recall that the Cauchy 
problem for the Navier-Stokes equations on M" is 

{ dtu -I- u ■ Vit — vAu -I- Vp = 0, 
div u = 0, 
u{0,x) = uo{x), 

where u : R x R" —>■ R" is the velocity vector field, ;/ > 0 is a viscosity parameter 
and p : R X R" —>■ R is the pressure. From now on, we consider only the cases 
n = 2,3, which are of special interest from the standpoint of classical mechanics 
and have been studied extensively over the past century. 

Among the first fundamental contributions to this problem, going back to 1934, 
were the works of Leray [20] and Hopf [16], establishing that ([HD admits global 
weak solutions 

(1.2) u e L°°([0 ,oo),L 2(K")) nL2([o,oo),R^(R”)) 

satisfying the energy inequality 

t 

(1.3) ||«(t,.)||i.+4 J ||Def(u(s,.))||i2d.< lluolli., 

0 
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where Def(u) = i(Vu+ (Vu)‘) is the so-called deformation tensor. Solutions of 
(EH) satisfying (USD and are said to be of Leray-Hopf type. 

When n = 2, it is known that Leray-Hopf solutions are, in fact, smooth and 
unique, but this is not known to hold for n = 3 (see [19], [7] for these and related 
results). 

In 1976, Heywood m considered (HD on a domain D C R”, under suitable 
boundary conditions, and showed that the uniqueness problem is related to the 
geometry of D. For instance, on 

(1.4) H = {(x, y, z) € : a; ^ 0, or < 1} 


there are nonunique solutions, whereas on bounded domains uniqueness always 
holds. 

In view of these developments, it is natural to consider the more general question 
of uniqueness of solutions of (HD on a general Riemannian manifold {M, g). 

In what follows, we are concerned only with the case when {M,g) is not flat. In 
this context, the system HD takes the form 


{ dtu + — vCu -\- Vp = 0, 

div M = 0, 
m(0,x) = uo(x), 

where V is the Levi-Civita connection, and C is an operator generalizing the Eu¬ 
clidean Laplacian (see below). When zz = 0 in the first equation of (11.51) , we obtain 
the well-known Euler equations 


( 1 . 6 ) 


dtu + VuU + Vp = 0, 
div w = 0, 
m(0,x) = uo{x), 


Throughout the literature, we encounter several natural choices for the operator 
C. The first one is the Bochner (or rough) Laplacian, given by As = V*V, where 
V* is the formal adjoint of V. Another possibility is to take C to be the Hodge 
(or damped) Laplacian (cf. |2]), given by Ah = dd* -|- d*d, where d and d* are 
the exterior derivative and its formal adjoint. Finally, one can also take (see m) 
Ah = 2Def*Def, where Def is the deformation tensor. 

The operators above, when acting on a divergence-free 1-form w, are related by 


(1.7) Ah{w) = Ab{w) — Ric(w, •) = Ah{w) — 2Ric(w, •), 


where Ric is the Ricci tensor and we identify 1-forms with vector fields via g. 

According to Ebin and Marsden m, when {M,g) is Einstein (i.e., Ric = Xg, 
for some constant A), the appropriate operator is C = Ad- Nevertheless, the other 
choices also appear in the literature. In particular, the system (11.511 with C = Ah 
(referred to in [6] as the modified Navier-Stokes equations), which we consider here, 
was studied in a number of works (e.g., i, m, 0, m, i)- 

The case of compact Riemannian manifolds {M,g), with or without boundary, 
has been studied by several authors (e.g., 0, m, m. [H])- For instance, in m 
it was proved that (HD is well-posed on Sobolev spaces, in any dimension. 

For non-compact Riemannian manifolds, however, the situation is quite differ¬ 
ent. In sharp contrast with the classical results in mentioned above, Chan and 
Czubak [6] proved in 2010 that there are nonunique Leray-Hopf solutions of (11.51) 
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on the two-dimensional hyperbolic space H^(—a^). In [6], the authors use C = Ad 
(since EI^(—a^) is Einstein), and also prove a nonuniqueness result where C = Ah 
for more general negatively curved two-dimensional manifolds. 

The question was raised in as to whether their methods could be extended to 
three-dimensional manifolds. 

In the case of EI^(—a^), the above question was answered negatively by Khesin 
and Misiolek [18] in 2012, who showed that the nonuniqueness phenomenon ob¬ 
served in [5] is essentially due to the Hodge-Kodaira decomposition, and that sim¬ 
ilar techniques cannot produce nonunique solutions on ]HI"(—a^), for n > 2. Some 
of their findings can be summarized in the following theorem, which incorporates 
an earlier result of Dodziuk |9] . 

Theorem 1.1. Given a Riemannian manifold {M,g), every harmonic vector 
field that belongs to the Sobolev space W^’'^(TM) defines a time-independent so¬ 
lution of the Euler equations (HD). Furthermore, there is an infinite-dimensional 
space of such vector fields on the Hyperbolic plane a^). On the other hand, 
for n > 3 any harmonic vector field on a^) vanishes identically. 

It is worth pointing out here that the space W^'^iTM) does not appear explicitly 
in [18] due to a typo in the statement of Proposition 1.2 of that reference. As one 
sees directly from its proof. Proposition 1.2 of [18] requires precisely this decay. 

The relationship between Theorem 11.11 and the Cauchy problem for the Navier- 
Stokes equations (1131) is as follows. Let u = be a time-independent solution of 
the Euler equations dm, for some function </>: M —)► R. If u decays as in Theorem 
11.11 and has finite dissipation (i.e., finite norm), one obtains infinitely many 
Leray-Hopf solutions u = f{t)v of (11.511 . for carefully chosen real-valued functions 
f{t). Because of the finite dissipation condition, one ends up having to require both 
and norms to be finite. 

This was the basic technique employed in [^ to obtain a nonuniqueness result 
in a^). In particular, one sees from Theorem 11.11 that the methods of [6] do 
not yield nonuniqueness of (II.5p in the case of a^). 

Nevertheless, we show that the approach introduced in [6] can be implemented 
on certain three-dimensional manifolds of non-constant negative curvature. 

As a matter of fact, the manifolds introduced by Anderson [1] in his counterexam¬ 
ple to the Dodziuk-Singer conjecture can be shown, after some delicate estimates, 
to carry nontrivial harmonic vector fields with sufficient decay: 

Theorem 1.2. The Anderson manifold {M^,g) has the property that the space of 
harmonic vector fields in W^’'^(TM) r\W^’'^(TM) is infinite-dimensional. 

We emphasize here that, in order to prove Theorem 11.21 it is not enough to 
repeat the construction of the harmonic forms as in [1]. We first need to produce 
a new, nontrivial supersolution to a certain PDE (system (12.15^ 1. This is done in 
Lemma 12.11 and provides stronger decay properties for the solutions of the given 
PDE. After that, we show that these decay properties are sharp ('Lemma 12.2p . in 
the sense that they are exactly what is needed for the harmonic forms to lie in the 
space 1T^’^(T*M) n W^’^{T*M). Finally, we carry out the Sobolev estimates in 
full detail. 

As an immediate corollary of Theorems 11.11 and 11.21 and the construction in 
[^ mentioned above, we obtain the following nonuniqueness result for the Navier- 
Stokes equations. 
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Theorem 1.3. The Cauchy problem (ll.5|l . with C = Ah, is ill-posed on 

in the sense that there are smooth initial conditions uq for which one has infinitely 

many smooth Leray-Hopf solutions of C31). 

One could state a version of Theorem ll.SI for the Euler equations (11.61) . However, 
our solutions do not conserve energy (i.e., the norm is not constant), and this 
is usually regarded as an essential feature of solutions of the Euler equations. 

The paper is organized as follows. In Section 2, besides reviewing the construc¬ 
tion of the manifold we prove the existence of a new supersolution that is 

more refined than the one used in [1] to create nontrivial harmonic 1-forms. 

The proof of Theorem 11.21 is completed in Section 3, where we carry out the 
estimates necessary to prove that the 1-forms constructed in Section 2 satisfy the 
hypothesis of Theorem 11.11 

For other results concerning the Navier-Stokes equations on a Riemannian man¬ 
ifold, see [6] and references therein. 

Acknowledgements. The author wishes to thank Prof. Gerard Misiolek for sug¬ 
gesting the problem and for many inspiring discussions, and Prof. Frederico Xavier 
for his interest and helpful comments on this work. This research was supported 
by the Richard and Peggy Notebaert Fellowship. 


2. The construction of { M ^, g ) 

Throughout this paper, whenever if —> M is a vector bundle equipped with a 
Riemannian metric, we denote by U’{E) the space of p-integrable sections of E and 
by IT®’P(if) the LP-based Sobolev space of sections of E, with smoothness index s. 
See [23] for more details. 

In this section, we recall Anderson’s construction of a complete Riemannian 
three-dimensional manifold {M^,g) and perform the modifications needed in our 
case. The latter are necessary in order to ensure that the constructed harmonic 
vector fields obey appropriate decay conditions at infinity. Our basic reference for 
this section is [1]. 

Consider the set 

(2.1) H = {{x,y) e C : < 1, a; > 0}. 


Fix a number a > 1 and equip H with the hyperbolic metric ds^ of sectional 
curvature —a^. In cartesian coordinates z = x iy, it takes the form 


( 2 . 2 ) 


dsi = 4 


dx^ -I- dy^ 

a2(l- |2|2)2- 


Consider a warped product H x f S^, where is the circle and / : i? —>• K is a 
smooth function to be defined. Let L be the line segment from (0, —1) to (0,1) and 
let Pa{z) be the distance, with respect to from the point z € H to L. More 
explicitly, we have (see (U p. 162]) 


(2.3) 


Paiz) = — arcsinh 
a 



Define / by 
(2.4) 


f{z) = sinho (pa( 2 ;)) . 
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The explicit formula for / will be used later on, but for now we only remark that 
from the above expression it follows that / extends smoothly to H U L, and that 
the corresponding warped metric extends to M = H x f S^, which is topologically 
just (for the metric, see [24l p. 12]). Using cartesian coordinates z = x + iy on 
H and the standard angular coordinate 0 on 5”^, the metric on M reads 

(2.5) = + 

To obtain a large space of harmonic 1-forms on M, we first decompose the Hodge 
Laplacian Am acting on 1-forms on M as 

(2.6) Am = doip + ip o d, 

where F = Vlog(/), ip denotes interior multiplication by F and V is the gradient 
of the Riemannian metric on M. Since / = 0 on L, the above formula only makes 
sense on H x f S^. However, if we require that solutions of the equation 

(2.7) Ah^(-.=)(w) - d(w(F)) + (dw)(F,.) = 0, 

satisfy a suitable Neumann condition at L, then they can be extended smoothly to 
M, and the extension will be harmonic. First, we require w to be invariant under 
isometric reflection through L. Second, we require that uj be invariant under the 
isometric action of the S'^-factor. This reduces the problem to solving equation 
(iTTll on H. 


Remark. From this point on, we specialize to the case a = 2 for simplicity. Other 
values of a > 1 would work as well, with minor modifications. 


Assume now that w = du. Let A = A]h[ 2 (_ 4 ). Then (12.7|) becomes 

(2.8) d{Au-du{F)) = 0, 


so that it is enough to solve 


(2.9) Au-du{F)=0 

on iJ. We would like to rewrite this equation on the strip 

(2.10) H = {(r, s) € C : 0 < s < 7r/2}, 

where it takes a simpler form. Thus, we introduce the biholomorphisms (/? : H 
and = -ip : F[ ^ fl given by 

- 1 

(2.11) ip{w) = pj{z) = log 

Note that ipi^) = R x {0}, and that the image of the full unit disk under ip is 

(2.12) O = {(r, s) : —7r/2 < s < 7r/2}. 

Letting v = uo (f, dSD) can be rewritten on H as 



d'^v d'^v 

Lniv) = ^ ^ 


^(log(/)) dv 

ds ds 


= 0 , 


(2.13) 
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This means that m is a solution of (12.9p if and only if u is a 
A long but straightforward computation using (j2.4p gives 

j, . 1 (1 + sin(s))^/^ - (1 - sin(s))^/^ 

o-rTT72-■ 

2 cos(s)^/^ 

Observe that / depends only on one variable. Each solution of (12.131) gives a 
harmonic 1-form on H. Consider the mixed boundary value problem 

(Lniv)=0, 

(2-15) J|^(^,o) = 0, 

[u(r,7r/2) = Voir), 

where vq : M —>■ M is a smooth function with uniformly bounded derivatives of 
all orders. It was shown in [T] that (12.151) has a unique solution v which can be 
extended by reflection to 

(2.16) O = {(r, s) : —-Kj^ < s < 7r/2}. 

This solution is smooth in Q and continuous up to the boundary 

(2.17) dQ = {(r, s) : s = 7r/2} U {(r, s) : s = —7r/2}. 

In particular, by construction, u = Uq on the boundary of 17. By the maximum 
principle applied to the operator Lq (see [H p. 35]), constructing appropriate 
supersolutions of (12.151) yields bounds on v and its first and second derivatives. 
More precisely, if is a positive supersolution of (I2.15P and v is chosen such that 
v{r,Tr/2) < u+(r,7r/2), then 

(2.18) jwj -I- \Vr\ + |Ws| + \Vrr\ + ^rs] + < const, v'^ 

holds pointwise. It should be remarked that (12.181) relies on estimates obtained 
from the method used to prove the existence of solutions to (|2.15p in [1], rather 
than a general argument. 

Let 

(2.19) ^■'■(r, s) = e“'*l’’lp(s), 

where p : [0,7r/2] —(0, -boo) is a smooth function and <5 > 0. We can see directly 
from (12.131) that v~^ is a positive supersolution of (12.151) (with initial condition 
U()’(r) = e“‘^l’'lp(7r/2)) if and only if 

(2.20) + (f) - “• 

.p'(0) = 0. 

In [T], the author uses p(s) = — cis^ -I- C 2 , for carefully chosen ci,C 2 and <5 > 0 
sufficiently small. In that case, for optimal ci and C 2 , the number <5 has to be less 
than ~ 0.78. However, our construction requires 6 to be at least 1, as shown 

by Lemma [2^ For this, we need the following 

Lemma 2.1. The function 

(2.21) u+(r, s) = e“l’'l(8s^ - -b 75), 
is a supersolution of (12.151) . 


where f = f o (p. 
solution of (12.131) . 

(2.14) 
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Proof. We must check that p{s) = 8s^ — 50s^ + 75 is positive on [0,7r/2] and that 
it satisfies (12.201) with (5=1. Positivity follows from the fact that the roots of p{s) 
are ±-^5/2, ±-\/T5/2, and p(0) > 0. The inequality in (12.201) becomes 

(2.22) (83"^ + 463^ - 25) + G(s)(32s^ - IOO 3 ) < 0, 


where 


(2.23) 


G(3) 


f / f^i ( \/l + sin(3) + ^l-sin(s) 
2 ^ 1/1 + sin(3) - y/1- sin(s) 



It is clear that G{s) > 0 and ^ 1 ( 3 ) = 323^ — IOO 3 < 0 on (0,7r/2). Only ^ 2 ( 3 ) = 
83 '^ + 403^ — 25 changes sign, so we break the inequality (12.221) in two cases: 


j 0 < 3 < V2/2 

\V2/2 < 3 <7r/2 

In the first case, P 2 {s) is negative, so (12.221) obviously holds. In the second case, we 
have 1/2 < sin(3) < 1, and this gives 


(2.24) 


G{s) > - ( 1 + 


2 cos(3) 


so that 


(2.25) 


(83^ + 463^ - 25) + G(s)(3233 


100s) < (83^^ + 4032-25) 


+ ( 103 ^ 


50s) + 


8s 3 - 25s 
cos(s) 


and this last expression is negative on the interval (•\/2/2,7r/2). In either case, 
(I2.22|) holds, thus proving the Lemma. □ 


We will only work with solutions v of (12.151) obeying (12.181) . Lemma [2+] implies the 
pointwise estimate 

(2.20) |n| + |n,.| + |ns| + |nrr| + |frs| + kssl < const. 

We finish this section with two Lemmas. The last one is the main ingredient needed 
in the next section to establish that du € W^’'^(T*M) D W^’‘^{T*M). 

Lemma 2.2. Let z € H and 'tp{z) = tpi{z) + i'tp 2 {z), where iPit' 4’2 o,re real-valued 
and ip is the function defined in (I2.11|) . The functions 

(2.27) /i(z) = (1 - |z| 2)|'!/;'(2)| and / 2 (z) = 

are bounded on H, provided that S is at least 1. If S < 1, then fi is not bounded. 


Proof. That fi is bounded follows directly from the triangle inequality. As for / 2 , 
note that |i/’^(z)| has only two singularities, located at ±i. We analyze the one at 
z = i., the other one being completely analogous. We have 


(x2 + (1 - y)^Y 
(1 + ’ 


(2.28) 


e 2'5 |V’i(^)I|.^'^2)|2 = const. 
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and switching to polar coordinates (r, 0), 


g-25|V.i(z)||^/(^)|2 

(2.29) 

(1 — 2r sin(0) -|- r^) 
-const. i^2r2cos(20 )+r' 

1 

« const. — -. 

(1 - r)2-2^ 

since 9 fs 7r/2, r « 1. □ 



const. 


(l-r)2^ 

(l + r)2(l -r)2 


Lemma 2.3. For u = v o ip, where v is a solution of (12.151) satisfying (I2.26|) . we 
have the pointwise estimates 


(2.30) \ux \ + l^yl < const, 
and 

(2.31) F “t" COnst. ^Ip |. 

Proof. From Lemma [2^ and (12.261) . we have 

(2.32) \ux\ + \uy\ < const. e~^'^^^\'ip'\ < const. 

For the second derivatives of u, using the Cauchy-Riemann equations and the bound 

\r{z)\ < w{z)\\ 

lUxxl = \'^rr{lpl)x F ‘^Vrs{lpl')x{lp2)x “t“ Vss{lp2)x F Vr{lpl) XX H“ '^s {.i’2) XX I 

< const. {{iplf F 2 \iPIiPI\ + {iflf + IV'Ll + I^Ll) 

(2.33) < const, {W?F\ip"\) 

< const. e~^'^ ^ \ip'\'^ 

< const. \ip'\. 


The same estimate holds for \uxy\ and \uyy\. □ 


3. Sobolev estimates 

In this section, we show that the and norms of du are finite. Recall 

that, up to a set of measure zero, M = H x f . Since is compact and the 
integrands will not depend on 9, we can reduce all computations to iJ. 

Denote by || ||g the norm with respect to the metric {gij) of M, so that from 
dlSl) we have 

\\dx\\g = \\dy\\g = {l-\z\^) and ||d6»||g = 1//, 

where z = x F iy as before. Then, using (12.301) . 

(3.1) J ll^wllg /= Ji^lF Ug)/ dx dy < const, j f dx dy. 

H H H 
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estimate. We show explicitly that the last integral in (jS.ip is finite. Expanding 
(|2.4p gives 


f{x,y) = - 


1 


2x 


+ 1 1 - 


Ax^ 


- 1/2 


1 


2x 


The first term above is bounded by 2 in absolute value, since 

2x 


+ 1 1 - 


Ax^ 


1/2 


(l-|zP)2 


(3.2) 


l-|zP 

As for the second term, note that 


+ 1 /1 ' 


Ax 


2x 


(3.3) 


1-UI 


+ 1 1 + 




(i-kIT 

1/2 


(1-|ZP)2 


< 


< 


> 1 . 


Ax 
^ Ax 


+1 


1/2 


+ 1 . 


Using polar coordinates. 


(3.4) 


f dx dy < const. + const. 


U2 fi 2r 


H 


I-7T/2 JO 


VA - 1-2 


drdO < oo. 


Therefore, we have du G L^{T*M). 


and estimates. From (13.41) and the fact that 

(3.5) lldull^ = (1 - \z\^nul + uir < const. (1 - 

we see immediately that du G L^{T*M). It remains to show that 

Vdu G L^{T*M 0 T*M) D L'^{T*M 0 T*M), 

where V denotes the Levi-Civita connection acting on 1-forms. First, we compute 
and find an upper bound for ||V(iu||g = (Vdit, Vdu). To do this, we break up the 
expression into three terms: 

II Vdullg = \\uxxdx (Si dx + Uxydx 0 dy -I- Uxydy (S dx + Uyydy 0 dy\\g 

+ 2(^Uxxdx 0 dx + Uxydx 0 dy + Uxydy (S dx + Uyydy 0 dy, UxVdx + UyVdy 
A- ||ua:Vdx -I- UyVdyWl 
and then estimate each one separately. 


Let us start from the last term containing only the lower order derivatives of u 
appearing in Vdu. Recall that, using Einstein’s summation convention, we have 

Vdx^ = -rf^ dx* 0 dx\ 

where are the Christoffel symbols, given in terms of the metric {gij) of M by 
the well-known formula 

rfc = \„ki 

2 ^ \dxj dxi dxi ) 
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From this we compute 


2x 2y 2y 

vdx = z - r-^dx 0 ax + - r~^dx 0 ay + --® 


1-kP 

—2x 


i-|zp 

dy 0 dy - (1 - \z\^Y f f^dO 0 d6», 


(3.6) 


Vdy = 


-2y 
l-|zP 
2 y 


dx 0 dx + 


2x 


jdx 0 dy + 


2 x 


l-|z|2 ^ l-|z|2 

dy 0 dy — (1 — \z\^)^ f fydO 0 d9. 

2 


dy 0 dx 


1-kP 

To obtain an upper bound for ||Ma;Vdx + UyVdyWg, we split 

||ua;Vdx + UyVdyWg = ||Ma:Vdx||g + ||My Vdyllg + 2{uo,Vdx, UyVdy) 
further into three parts: 

J X 


Vdxll^ = ui[l - (^4x^ + 4y^ + 4y^ + 4x^ + (1 - \z\^Yip 

VdyWl = ulil - |z|2)2 Uy^ + 4x2 + 4x2 + 4y2 + (1 - \z\Y^^ 


2 1 ’ 


2{ua:'Vdx, UyVdy) = 2uxUy{l — \zf) ( —4xy + 4xy + 4xy — 4xy + (1 — \zf) 


|2^2 fxfy 


p 


= 2UxUy{l - |z|2) 
Adding all three terms gives us 


2\4 fxfy 


p 


llu^Vdx + MyVdy||2 = 8|z|2(u2 + - |z|2)2 


(3.7) 


2\4 


(1-kn 

p 




Now, recall that u satisfies (12.91) . that is, Au = du{F), where F = Vlog(/). In 
coordinates, this becomes 

(3.8) Au = dM(V10g(/)) = —Uxfx + ^Uyfy = - J - —{Uxfx+Uyfy), 

where (y®^ ) are the components of the inverse of y, so that 

/ 


(3.9) 


'^xfx T ^yfy — 


-Au. 


^vJy- ( i _| 2 | 2 ) 2 - 
Substituting this back into (13.71) gives 

(3.10) lla.Vdx + a,Vdy||2 = 8|z|2(a2 + ul){l - |z|2)2 + (Aa)2. 
Thus, it remains to estimate |Am|. Using (12.311) and Lemma 12.21 we get 

“ In '^xx T y ‘^yy\ 

(3.11) < const. (1 — |z|2 )2|')^'(2;)| 

< const. (1 — |z|2 ). 
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Going back to (13.1011 and using ()2..30|1 . we get 
(3.12) Wux^dx + UyVdyllg < const. (1 — 

This takes care of the term containing only first order derivatives of u in (Vdu, Vdu). 
For the term containing only second order derivatives, using Lemmaand (12.311) 
gives 

\\uxxdx 0 dx + Uxydx 0 dj/ + Uxydy 0 dx + Uyydy 0 dy\\^ = 

= {ulx + 2ul^ + ul^){l-\z\^r 

< const. |i/’'(2;)p(l — 

< const. (1 — 1-2^)^, 

which is the same bound as before. Finally, for the mixed term that contains first 
and second derivatives of u in (Vdu, Vdu), we again use (12.3011 . (12.3111 and Lemma 

[O 

^(uxxdx 0 dx + Uxydx 0 dy + Uxydy 0 dx + Uyydy 0 dy, UxVdx + UyVdy 
<2(1— \Ux \ ^2|Uxx| + 4|Uxy I + 2|Myy 1^ + I ^2|u3,2,| + 4|u3,y| + 2|Uyy|^ 

< const. (1 - |2p)^ \Ux\W{z)\ + \Uy\\lp\z)\ 

< const. (1 — |2|^)^|'!/’^(2)| 

< const. (1 — |2p)^. 

Combining all the above estimates and squaring leads to 
(3.14) llVdullg < const. (1 — \z\^Y, 

so that 


(3.15) 


< const. [ f dx dy 


H H 

and this integral is finite (see (13.41) 1. This shows that du £ 1F^’^(T*M). 
Moreover, one sees directly from the above estimate that 

dx dy 


(3.16) 


H 


{l-\z\2)2 


< 00 , 


so that finite dissipation, i.e., du G is also satisfied. Combined with the 

estimate, we have du G M) D M). 
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